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NUMERICAL SOLUTIONS OF FREE-MOLECULE FLOW 


IN CONVERGING AND DIVERGING 
TUBES AND SLOTS 

by Edward A. Richley and Thaine W. Reynolds 
Lewis Research Center 


SUMMARY 

The integral equations governing flux distributions for free-molecule flow 
through converging or diverging tubes and slots are developed under assumptions 
of uniform entering flux and diffuse particle reflection from the walls. The 
equations are solved on an IBM 7094 computer by applying an iterative method of 
solution. Specific cases studied in the analysis are tube length to radius 
ratios of l/2, 1, 2, 4 , 8, and 16; slot length to width ratios of l/4, l/2, 1, 
2, 4, and 8; and, for each ratio, up to 14 variations of wall half-angle rang- 
ing from 75° to -60°. Wall and exit-aperture flux distributions and total and 
direct transmission probabilities are presented as functions of these param- 
eters. Transmission probabilities for a wall half-angle of zero are compared 
with results of other investigations and show excellent agreement. Results are 
presented in tabular and graphic form, and reasonably accurate results for 
parametric values within the specified ranges may be determined by interpola- 
tion. 


INTRODUCTION 

Gas flow characteristics under conditions of free-molecule flow have been 
of interest since the early work of Knudsen and Clausing and such study has 
gained impetus in recent times primarily because of the rapid expansion of 
space technology. Applications of solutions of free-molecule flow include a 
wide range of problems, for example, the design of electric rocket thrustors, 
vacuum gages, and vacuum facilities (refs. 1 and 2). The subject of this 
report is free-molecule flow through converging or diverging tubes and con- 
verging or diverging two-dimensional slots. 

Mathematical formulation of such free-molecule "internal” flow problems 
was first accomplished by Clausing for right-circular tubes. Later investi- 
gators rederived Clausing f s equation by different techniques and also formu- 
lated the problem for other models (e. g. , ref. 3). 

The problem discussed herein is to determine the flux distribution along 



the walls of the tube or the slot for the case of random, free-molecular gas 
flow into the inlet. 

Particle reflection from the walls is considered herein as in most other 
investigations, to be diffuse. This is an appropriate assumption from a phys- 
ical viewpoint, if the particles are adsorbed and then evaporated from the 
wall. 


Of course, the presence of an adsorbed layer of particles on the surface 
implies the possibility of surface diffusion. Evaluation of the effect of 
surface diffusion on flow requires knowledge of the diffusion coefficient and 
equilibrium adsorption data for the particular system considered. For example, 
it is shown in reference 4 that, for silver atoms effusing through molybdenum 
and nickel orifices, the total flow may be 30 percent greater than that cal- 
culated for vapor flow alone. For the effect to be of this magnitude, it was 
required that the length-radius ratio of the orifice be less than 1, while the 
radius was small (R < 0.5 mm). In order to keep the analysis herein independ- 
ent of the gas or tube material, surface diffusion effects were not considered. 
A more detailed discussion of gas-surface phenomena is given in references 5 
and 6. 


Once the wall flux distribution has been determined, the remaining flow 
characteristics are readily calculated. The general expression for the wall 
flux distribution is an Integral equation, and solution in closed form is usu- 
ally not possible without an additional assumption. An assumption that has 
been used for the right- circular- tube problem, for example, is that particle 
flux on the tube wall varies linearly with distance from the entrance (ref. 7). 
As mentioned in reference 7, however, the integral equation is amenable to 
direct integration on a digital computer, and no assumption regarding the wall 
flux distribution is required. An iterative method of solution, on a digital 
computer, is described in reference 8 for the problem of the right-circular 
tube and in reference 9 for diverging tubes with wall half-angles ranging from 
0° to 22.5°. As discussed In reference 9, solutions obtained with the assump- 
tion of a linear wall flux variation for the right- circular tube are in approx- 
imate agreement with numerical solutions. This assumption, however, is no 
longer applicable for a tapered tube. 

An iterative numerical method of solution that is similar to those re- 
ported in references 8 to 10 is applied herein to cover a much wider range of 
geometric variations than were available previously. Wall and exit-aperture 
flux distributions and total and direct transmission probabilities are pre- 
sented in both tabular and graphic form. The computer program used to obtain 
these results is given In appendix C by Carl D. Bogart. 


ANALYTIC RELATIONS 

Configurations and Basic Flux Equations 

The configurations investigated herein are shown in sketch (a) and fig- 
ure 1. 
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Tubes 



Slots 

(a) 


Particle flow into the tube or slot is assumed to be uniform over the in- 
let and to have random distribution of direction. Explicit knowledge of the 
speed distribution function is not a requirement for this problem (see ref. 7). 
This implies that solutions are independent of the wall or particle tempera- 
ture. Downstream of the exit, vacuum conditions are assumed so that there is 
no return flow from this region. Free-molecule, or Knudsen, flow is assumed 
in the region internal to the tube or slot; thus, particle-particle collisions 
are negligible and only particle-wall collisions need be considered. Particle 
reflection from the walls is assumed diffuse, that is, assumed to follow the 
cosine law. With these assumptions, the particle arrival rate at any point on 
the wall may be determined numerically on a digital computer. 

The basic relation from which the particle flow behavior may be derived 
is (ref. 7) 


dn b 


n a cos e a b cos e^a 



*Aa 


(i) 


where 
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drift flux (particles per unit area per unit time) arriving at differential 
area dAft from dA a 


n a flux leaving differential area dA a 


^ab 

distance from 

a to b 



^ab 

angle between 

^ab 

and 

the normal to 

<3A a 

^ba 

angle between 

^ab 

and 

the normal to 

dA-b 


The geometric relations expressed in equation (l) are shown in sketch (b) 
and all symbols are defined in appendix A. 




(c) 


The total arrival rate nft at dAft 
is obtained by integrating equation (l) 
over all the source area A a that con- 
tributes to the flux at dAft. 

If the configuration being exam- 
ined is two-dimensional; that iS; of in- 
finite extent in the third dimension; 
such as shown in sketch (c); the par- 
ticle flux relations can be developed 
from the line-source relation 


dn b = 


n a cos e ab cos e ba 


27 




( 2 ) 


ab 


where 


dnft flux (particles per unit area per 
unit time) arriving at the dif- 
ferential area dAft from line 
element dy a 

n„ flux leaving line source dy n 
a a* 

Z a ft distance from b to line element 


Equation (2) is, of course; obtained from equation (l) by the partial integra- 
tion of dA a over the line element from z equals negative infinity to posi- 
tive infinity. The total arrival rate nft at dAft is obtained by integrating 
equation (2) over all line elements that contribute to the flux at dAft. 

The following consistent subscript notation is used in all subsequent de- 
velopments: 


Subscript 

Surface or plane 

i 

2 , 3 
4 

Inlet plane 
Wall surface 
Exit plane 
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Tube Configuration 


The flow relations for tubes are given in the following paragraphs. Some 
identities helpful in reducing the general relations are given in appendix B. 

Tube wall flux . - The flux ng(xg) on the tube wall at any point (see 
fig. 1 (a)) consists of particles arriving directly from the inlet plus par- 
ticles arriving from the remainder of the wall. When the general flux equa- 
tion (l) is employed, the flux at point 2 becomes 


n %(. x z) 



n l 


COS 0^2 COS 0oq 


Til 


2 

12 


dAq 



n 3 (x 3 ) 


cos 0 2 3 cos 0^2 


Til 


2 

23 


dA 3 


(3) 


Since n^ is constant, the first term of equation (3) may be integrated 
directly. Expressed in dimensionless form, it becomes 


1 

2 r 2 


J x| sec 2 p + 3xo tan p + 2 _ _ 

< / x Xg cos p - r 2 sin 3> 

yXg sec p 4- 4r 2 J 


(4) 


Barred quantities are nondimensional variables expressed as ratios to the inlet 
radius Rq. 


Since ^^ 3 ) is variable over the area A 3 and is also equal to the un- 
known arrival rate n 2 (x 2 ) that is being sought, the second term of equa- 
tion ( 3 ) cannot be completely integrated directly and must be finally deter- 
mined by some other means. A numerical method that employs an iterative tech- 
nique is used herein. First, however, since the flux n 3 (x 3 ) is independent of 
the angle a 3 (see fig. l(a) and appendix B) , integration with respect to a 3 
can be carried out to yield 


L 


n 3( x 3) S 1 - ( x 3 - x 2 ) sec p 


(x^ - Xo ) 2 sec 2 p + 6 r. 


2 I 3 


[(X 3 - X 2 ) 2 sec 2 p + 4r 2 r 3 ] 


3/2 t dx 3 


(5) 


The complete expression for n 2 (x 2 ) in nondimensional form is 


1 [ X 1 sec2 P + 3X2 tan P + 2 

n l 2r« 


/—2 2 

|rx 2 sec p + 4r g 


- xg cos P-Tg sin pi 


/ n 3^ x 3^ cos 3 L /- - x 


(x 3 - X 2 ) 2 sec 2 3 + 6r 2 r 3 


[(x 3 - Xg) 2 =— 2 


sec‘p + 4r 2 r 3 ] 


3 72 




( 6 ) 
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In this equation, as in subsequent ones, the half-angle 0 is positive for a 
divergent configuration and negative for a convergent configuration. 


Tube exit-plane flux . - The flux at a point in the exit plane (see fig. 
1(b)) consists of particles arriving directly from the inlet, plus particles 
arriving from the wall: 


nq(r 4 ) 



COS 614 . COS 041 

n l ^ 


nl 


14 



cos 6 cos e 42 



dA 2 


( 7 ) 


Integration of equation (7) expressed in nondimens ional form yields 



The first term of equation (8) represents the flux of particles arriving di- 
rectly from the inlet opening. The second term is the flux of particles that 
have experienced one or more collisions with the wall before arriving at a 
location r^ in the exit plane. 

Tube transmission probability . - The fraction of particles incident upon 
the tube inlet plane that passes through the tube without colliding with the 
wall, that is, the direct transmission probability Pq, is obtained by inte- 
grating the direct flux portion of equation (8) over the exit area of the tube. 
The resulting expression for P^ is 


P 


d 


1 

2 


_2 _2 

1 + L + R L 



( 9 ) 


The total transmission probability P^, or the fraction of particles in- 
cident upon the tube inlet area that eventually pass out through the downstream 
end of the tube, is determined by integrating the total flux nq over the exit 
area and dividing by the total inlet flux: 


Pt 



dr. 


( 10 ) 
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II 


Slot Configuration 

The flow relations for slots are given in the following paragraphs. Some 
identities helpful in reducing the general relations are given in appendix B. 

Slot wall flux , - The basic line source relation (eq. (2)) is used to ob- 
tain the arrival flux at the wall due to direct flow from the inlet and flow 
from other elements of the opposite wall such as shown in figure l(c): 


n 2 (x 2 ) = 


/ 

* / W n 


n l 


COS 0]_ 2 COS 0 2 ]_ 


21 


&Y-\ 


12 


/ 


n 3 (x 3 ) 


cos 0 2 3 cos 0 32 


2Z 


a# 


( 11 ) 


23 


It will be noted that in this configuration no flux arrives at point 2 from any 
other point on that same wall, since all such points are located at angles of 
90° from point 2. After the proper geometric relations are substituted into 
equation (ll) (see appendix B) and the indicated integrations are performed, 
the following nondimensional expression for the flux on the wall is obtained: 


ng ^ X2 ' > = i <( 1 - cos 0 
n ± 2 % 


+ 1 


/ 


x 2 sec^jB + tan (3 
L^x| sec 2 p + 2 x 2 tan p + 1_ 


n 5( x 5) 

n l 


cos p- 


^ Q 

1 + 2 (x 2 + X 3 ) tan p + 4 x 2 X 3 tan^p 
[( x 2 “ x 3 ) 2 + C 1 + x 2 tan p + x 3 tan p) 2 ] 


3721 


dx. 


( 12 ) 


The first term in brackets in equation (12) is the arrival rate at a point 
on the wall due to direct flux from the inlet. The second term in brackets 
represents the flux arriving from the opposite wall. Barred quantities are 
nondimensional variables expressed as ratios to the inlet slot width. 


Slot exit-plane flux . - The flux at a point y 4 in the exit plane (fig. 
1 (d)) consists of contributions from the open end directly and from the two 
walls. The equation expressing these terms is (see appendix C) 


_ 1 f 1 
n i ' 2 IV^T 


= *-= 
(1 - 2y 4 ) 2 


1 + ) 

+ (1 + 2y 4 ) 2 J 




W L + ^4 Wl - 2y 4 

[4(L - x 2 ) 2 + (1 + 2y 4 + 2xg tan j3) 2 ] 3 / 2 [^(L - Xg) 2 + (l - 2y 4 + tan p) 2 ] 3 ^ 2 


■} 


^2 


(13) 
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The -first term in brackets in equation (13) represents particle flux at a 
point y^ in the exit plane coming directly from the inlet without a wall 
collision. The second term in brackets represents flux of particles at 
point y^ that experience one or more wall collisions before arriving at 
that point. 


Slot transmission probability . - The fraction of particles that are in- 
cident upon the slot inlet and pass through the slot without colliding with 
the wall is P d and is obtained by integration of the first term of equa- 
tion (13) over the exit area: 

P d = |[VlL 2 + (1 + W l )2'- V^L 2 + (1 - W L ) 2 ] (14) 


The total transmission probability P^ is the integral of both terms of 
equation (13) over the exit area: 


p t - 



n 4 (y 4 ) 

n l 



(15) 


Additional Relations 


Wall fluxes. - As previously mentioned, the tube and slot wall flux rela- 
tions, equations (6) and (12), may be used for either converging or diverging 
configurations by applying the appropriate sign to (3, the wall half-angle. A 
relation between the wall flux values for flow through any particular configu- 
ration (sketch (d)) in the forward or reverse directions may also be developed 
from consideration of equilibrium requirements. 



(d) 


Consider the flow through the open- 
ing from each chamber (sketch (d)) in- 
dependently with the assumption that the 
other chamber is at zero pressure. The 
flux at some point on the wall at a dis- 
tance x from the inlet A d may be 
calculated first for flow from the left 
chamber and then for flow from the right 
chamber. Let these two wall fluxes be 
designated ( n x) d4 and ( n x) 4 p.> respec- 
tively. Since these flows have been 
assumed to occur under free-molecule 
conditions, with no interaction of par- 
ticles, the flows may be superimposed 
without affecting the separate solu- 
tions. The total arrival rate at the 


wall at point x is then (n x ) d4 + ( n x) 4 i* 
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If the gas in the total enclosure is considered to he in a steady-state, 
the arrival rates nq and n^, and indeed the arrival rates anywhere within 
the enclosure, must all he equal. Thus, it can he stated that the wall arrival 
rate at point x under these conditions is 


or 


n l = ( n x)l4 + (Oxhl 

(n x )41 _ (n x )l4 
n l “ “ n l 


(16) 


(17) 


This relation, which does not depend on the shape of the opening, or passage- 
way, provides a convenient tool for evaluation of reverse flow through 'a par- 
ticular configuration, (e. g. , tube or slot) once the flow in the forward direc- 
tion has heen'^calculated from equation (6) or (12). 

Transmission prohahility . - Under steady-state conditions in the enclosure 
in sketch (d), .there must he _ no net flow through the tube, or slot, connecting 
the two chamber^. From this steady-state requirement, a relation between the 
transmission probabilities through the same configuration in the forward and 
reverse directions can be derived. Again, consider the flows from the two 
chambers to be occurring independently of each other. The particle flow from 
chamber 1 to chamber 4 is nqA^Pqp)-;- ; the particle flow from chamber 4 to 
chamber 1 is npA 4 i 5 \P 4 p)q, where (PqqJt an & (^4l)t are transmission 

probabilities for the respective flows. When the flows are superimposed, the 
respective total transmissions are not affected. Under steady-state conditions 
then, since there must x be no net flow, 

n l A l( P 14)t = n 4 A 4 ( P 41 ) t (18) 

But, also under steady-state conditions nq - n^, so that equation (18) becomes 

( P 14 ) t = n ( p 4l)t (19) 


As long as free-molecule flow conditions exist, the values of transmission 
probability, as defined herein, do not depend on flow rate. Equation (19) is a 
general relation, then, dependent only on configuration parameters. 

Consideration of the direct flux relations (eqs. (9) and (14)) shows that a 
similar relation holds for the direct transmission probabilities as well: 

( P 14 ) d. = ^ ( p 4l)d (20) 


If, for example, the configuration variables for the tube are introduced 
into equation (19) or (20), the equation may be written as 
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( 21 ) 


where 


(P 14^ L B 

S ’ 13 



L_ 

L_ = Rj 

r l i + J±. m 
R 1 


( 22 ) 


(see fig. 1(b)). 

For the slot configuration equation (19) or (20) becomes 



where 


L_ 

L _ W 1 

W L 1 + 2 — M 
W n 


(23) 


(24) 


(see fig. l(d)). 

The transmission probabilities presented in the next section were explic- 
itly determined for various configuration parameters and flow directions. By 
application of equation (2l) or (23), values of the reverse-flow transmission 
probabilities may be readily determined. 


METHOD OF SOLUTION 

Numerical solutions for the preceding equations were obtained on an 
IBM 7094 computer. The FORTRAN program is given in appendix C. 

Solution of the wall flux equations (eqs. (6) and (12)) was accomplished 
by application of an iterative procedure similar to that described in detail in 
reference 10. The procedure, as applied to equation (6), for example, consists 
in supplying an initial guess of the unknown function n3 ( X3 ) in the integral 
of equation (6). The equation is numerically integrated by using a combination 
of Simpson* s rule and trapezoidal integration. New values of the function are 
thus produced that are used in turn, in the next iteration. Iteration proceeds 
in this manner until convergence is reached, that is, until the maximum change 
in any one of the final pointwise values, compared with its previous value, is 
less than some preassigned limit (for details see appendix C). For the solu- 
tions presented herein, the largest value used for this limit was 0.02 percent. 
Similar to the findings of reference 10, trial solutions in which widely dif- 
fering initial guesses were employed always resulted in convergence to practi- 
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cally the same final answers; however, the rate of convergence differed. The 
simple initial guess of ^( 3 ^ 3 ) = 0 generally yielded the most rapid rate of 
convergence. 

After the wall flux relations are solved, the exit-plane equations 
(eqs. (7) and (13)) may he solved by direct numerical integration. Similarly, 
the total transmission probability equations (eqs. (10) and (15)) are amenable 
to direct numerical integration by use of the exit-plane flux values. 

The accuracy of the results, of course, depends on the size of the in- 
crement used in the numerical integration of the various equations and, in 
general, improves as the increment size is made smaller. Smaller increments, 
however, increase machine computation time, particularly with respect to the 
wall flux equations. Consideration of both these factors led to the selection 
of varying increment sizes such that the total number of calculated points was 
approximately the same for a.1 1 configurations. 


RESULTS AND DISCUSSION 
Transmission Probabilities 

The transmission probabilities calculated by equations (9), (10), (14), 
and (15) are presented in tables I and II. A comparison of the calculated 
total transmission probability for the straight tube and the parallel- walled 
slot (i. e. , for 0=0) with that of previous investigations is shown in fig- 
ure 2. The results of the present calculation’ by the iterative numerical 
technique are in agreement with those of reference 3. The transmission proba- 
bilities for 0=0 obtained in reference 3 were determined analytically under 
the assumption of a linear variation in wall flux. The numerical procedure 
used herein, as well as that used in reference 8, shows that this assumption 
is quite satisfactory for this case; however, it is not applicable when the 
wall half-angle departs from zero. 

Total transmission probabilities for the convergent and the divergent 
tubes are shown in figure 3(a) for length to inlet-radius ratios up to 16. For 
the divergent tubes, the transmission probability appears to become asymptotic, 
especially for half -angles greater than about 20°, where the total transmission 
probability levels off at a length to inlet-radius ratio of about 10. The re- 
sults that may be compared are in agreement with graphical results given in 
reference 9, which covers positive wall half-angles from 0° to 22.5°. Total 
transmission probabilities for the convergent and the divergent slots are 
shown in figure 3(b) for length to initial- width ratios up to 8. The behavior 
is qualitatively similar to that of the tubes. 

The direct transmission probabilities P^ for the tubes and slots are 
shown in table II and figure 4. These probabilities rapidly approach limiting 
values as the length- radius ratios, or length- width ratios, approach infinity. 
The values sin^0 for the tube and sin 0 for the slot are the limits of 
equations (9) and (14), respectively. Physically, as the ratio of length to 
initial width (or radius) increases for a fixed wall half-angle, the percentage 
of approaching particles that can escape directly through the exit at angles 
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greater than the wall half-angle (3 tends to decrease as shown in sketch (e). 


The fraction of in- 
cident particles that 
experience one or more 
wall collisions before 
passing out of the tube, 
or slot, is simply equal 
to P -5 - P< 3 _. This frac- 
tion may be determined 
from the values given in 
(e) tables I and II and is 

compared for the case 
where JB = 0 with total 

and direct transmission probability values in figure 5. The P^ - P^ curve 
passes through a maximum for this case as the length to inlet-radius ratio or 
length to inlet-width ratio increases. This behavior of the Pt - Pd curve is 
to be expected if the total and the direct transmission probabilities each 
approach zero with increasing length- radius or length-width ratios as they do 
for zero or negative wall half-angles. For the divergent configurations, how- 
ever, the total and the direct transmission probabilities do not approach zero 
(figs. 3 and 4) and a maximum in the Pt - Pd curve does not generally occur. 

Equations (21) to (24) may be used in connection with the data of tables I 
and II to obtain reverse-flow transmission probabilities for the various con- 
figurations. For example, consider a slot of length to inlet-width ratio of 
l/2 and a wall half-angle equal to 45°. The reverse-flow- configuration param- 
eters are a wall half-angle of 45° and from equation (24), a length to exit- 
width ratio of 1/4. From table 1(b) 

Pt = ( p 14) L „ = °- 986 
w l 

For flow in the opposite direction, equation (23) is used to calculate P-t: 

Pt=( p 4l)L =0.493 

This value, calculated from equilibrium requirements, is identical with the 
value given in table l(b) that was calculated from the theoretical flow equa- 
tions. 



Wall Flux Distributions 

Tubes . - Values of the incident flux or arrival rate at the wall of the 
convergent or divergent tube are given in table Ill(a). Several illustrative 
plots of data from this table are shown in figure 6 . Wall flux varies nearly 
linearly along the tube length for (B = 0 , is concave downward for the con- 
verging tube, and is concave upward for the diverging tube. 
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Slots , - Wall flux values for the convergent or divergent two-dimensional 
slots are presented in table 111(b) , and several illustrative plots are shown 
in figure 7. The wall flux distributions for the slots follow qualitatively 
the same behavior as those for the tubes. 

This behavior of the tube and slot wall flux distributions with varying 
wall half-angle is due to the varying magnitude of contributions received 
directly from the inlet and contributions from particles reflected from the 
walls. 

The qualitative behavior of these components of the flux curves is shown 
in sketch (f) for a length to inlet-width ratio of 1 and wall half-angles of 
± 20 °. 



(f) 


The wall fluxes presented herein for both the tube and the slot configu- 
rations may be used in equation (17) to obtain flux values for reverse flow 
through each of the configurations. These correlations are not obvious in the 
plots of figure 7, since the dimensionless parameter length-radius or length- 
width ratio given thereon is expressed in terms of the flow entrance aperture. 
Equations (22) or (24) should be used along with equation (17) to determine 
the appropriate reverse-flow configurations. 

The accuracy of the flow equation calculations may be checked by using 
some results given in table III. For example, length to inlet-radius ratios 
and wall half-angles of 1 and 45°, and l/2 and -45 , respectively, result in 
similar tube configurations. The calculated values from table III (a) (rounded 
off to three places) are reproduced and summed for comparison in the following 
table: 
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Axial distance. 

Length to inlet-radius ratio, 

Sum of flux 

x/L or 1 - x/L 

L/Ri 

ratios 


1 

1/2 



Wall half-angle, 



0, 

deg 



45 

-45 



Flux 

ratio. 



n 2/ n l 

1 


0 

0.166 

0.834 

1.000 

.1 

.135 

.865 

1.000 

.2 

.111 

.889 

1.000 

. 3 

.093 

.907 

1.000 

.4 

.078 

. 922 

1.000 

.5 

.066 

. 934 

i. ooo 

.6 

.056 

. 944 

1.000 

. 7 

.048 

. 952 

1. 000 

.8 

.042 

. 958 

1.000 

. 9 

.036 

. 964 

1. 000 

1.0 

.031 

.969 

1. 000 


The equilibrium requirement of equation (17), that the sum should equal 
1. 000 at each point, is thus satisfied. 


Exit-Plane Flux Distributions 

The integrals in the exit-plane flux distribution equations (eqs. (8) and 
(15)) for both the tubes and the slots contain functions that become somewhat 
difficult to evaluate accurately by numerical methods at the exit plane near 
the wall. The problem may be overcome to a great extent by use of very small 
increments; however, practical considerations arise in deciding on an appro- 
priate increment size. Results were determined to be fairly accurate over at 
least 95 percent of the exit opening. There is a possible error over the re- 
maining 5 percent, and this possibility is indicated in the figures by the 
dashed portions of the curves near the walls. 

Tubes . - Flux values across the exit plane at different radial distances 
from the tube axis are presented in table IV(a) for the convergent and the 
divergent tubes. Flux distributions across the exit plane of the cylindrical 
tubes of various lengths are shown in figure 8. The distributions tend to 
become flat with increasing length to radius ratio. 

Flux distributions across the exit plane of tubes of various wall half- 
angle are shown in figure 9 for length to inlet-radius ratios of 0.5, 2, 
and 16. As would be expected from the transmission probability results, the 
magnitude of the curves decreases with increasing length to radius ratio. 

Slots . - Flux values across the exit plane at different lateral distances 
from the centerline of the slot are given in table IV(b) for the convergent 
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and the divergent slots. Flux distributions at the exit plane of the parallel- 
walled slots of various lengths are shown in figure 10, and several illustra- 
tive plots of the variations of exit-plane flux distribution with wall half- 
angle are shown in figure 11. As is true of the wall flux curves , the general 
behavior of the exit-plane curves is qualitatively similar for slots and tubes. 


CONCLUDING REMARKS 

The object of this report was to determine the flow characteristics of 
converging or diverging tubes and slots under conditions of free-molecule 
flow. An iterative numerical method of solution of the integral equations 
that describe the flow was employed. With this method, no a priori knowledge 
of the flux distribution along the wall was required, as it would be if closed- 
form analytic solutions were sought. Wall flux distributions, exit-plane flux 
distributions, and total and direct transmission probabilities were determined 
for a wide range of configuration length to radius (or, width) ratio and wall 
half-angle. Relations were also developed that may be used to determine the 
"reverse-flow” wall flux distribution and the transmission probability of a 
given configuration directly from the computed "forward-flow" results. 

It was found that values of transmission probability for the configura- 
tion having a wall half-angle of zero were in agreement with values determined 
by other investigators. Some values presented for other configurations (i.e. , 
wall half-angle not equal to zero) were found to be in agreement with graphic 
results of another investigation. 

Wall flux distributions for the tubes and the slots were qualitatively 
similar to each other and were found to vary nearly linearly with distance for 
the zero wall half- angle configuration. For the divergent configurations the 
wall flux varied more sharply with distance near the entrance, while for con- 
vergent configurations it varied more sharply near the exit. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 11, 1964 
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APPENDIX A 


SXMBOLS 

A area 

L length of configuration, measured normal to entrance plane 
length of slot vail, measured along vail 
Z length of line, see sketch (b) 

M tangent of angle (B of tube or slot 

n flux or arrival rate, number/ (unit area) (unit time) 

Pp direct transmission probability, fraction of entering particles that exit 
downstream without a vail collision 

Pp total transmission probability, fraction of entering particles that exit 
downstream 

Rp inlet radius of axi symmetric configuration 

Rl exit radius of axisymmetric configuration 

r radial distance 

s length of surface normal from point on surface to centerline of configu- 
ration 

t distance, for example, t-y, distance from point i to end of s j that is 
on centerline 

Wp width of exit, two-dimensional configuration 

Wp width of inlet, two-dimensional configuration 

x axial distance 

y lateral distance 

a position angle in cylindrical coordinate system 

p wall half-angle of tube or slot, positive for diverging and negative for 
converging configuration 

6 angle between surface normal and line Z, see sketch (b) 
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Subscripts: 

a,b general points 

1 inlet plane 

2 , 3 wall 

4 exit plane 

Superscript: 

( ) dimensionless quantity for tube configuration when divided by 
for slot configuration when divided by W^_ 


1 anc ^ 
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APPENDIX B 


IDENTITIES FOR INTEGRATION OF EQUATIONS 


The following identities for the tube configuration may be noted from 
figures l(a) and (b): 


dA]_ = r^_ dr^_ dcx^ 


„ x 2 

cos e 12 = — - 

z 12 


l \ 2 + s| 

cos 021 = —== 

2 2Z 12 s 2 


- t 


12 


t 2 g = r| + (x 2 + r 2 tan p)‘ 


l lz = r ! + r 2 - 2r l r 2 cos a l + *2 

(IA3 = r3 1x3 dcc.3 sec |B 


,2 ^ o 2 + 2 

„„ n _ z 23 + s 3 - ^23 
cos e 32 — — — 

^23 s 3 


2 2 2 

rn „ P _ Z 23 + s 2 “ t 32 
cos 0 23 — — — 

^23 2 


t| 2 = r| + (x 3 - x 2 - r 2 tan p) 2 

t 23 = r 2 + ( x 3 - x 2 + r 3 tan p) 2 

s 3 = r 3 sec 0 

s 2 = r 2 sec p 

r 3 = Rl + x 3 tan P 
r 2 = Ri + x 2 tan p 

Z 23 = r 2 + r 3 " 2r 2 r 3 cos a 3 + ( x 3 ~ x 2) 2 
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cos 0j_4 = cos 0/n = 


41 


Z 14 

dA^ = rg dag dxg sec p 


,2 ^ ^2 


^ 14 = r i + r 4 - 2r l r 4 cos(a 1 - a 4 ) + if 


z 2 4 = r | + r 4 " 2r 2 r 4 cos a 2 + ( L “ x 2^ 2 

L - x 2 


cos 0 42 = 


z 24 


7 2 4. o 2 4-2 

cos 0o >. = _ 2 42 


^24 = 


2Z 24 s 2 


t 4 g - (L “ X 2 ~ r 2 ' tan 0) 2 + 


2 , 2 


The following identities for the slot configuration may he noted from 
figures l(c) and (d): 

d-S? = <1x3 sec 0 
cos 0 12 = 

z 12 


cos 0 


21 = 


7 2 , 2 ,2 

z 12 + s 2 ~ ^12 

2 z 12 s 2 


s 2 = ^2 sec P 


^12 = yf + (x 2 + y 2 tan 0)' 

z i2 = x l + (y2 “ yi) 2 


cos 0 23 = Z 25 + S 2 - t2 2 
^ 2Z 23 s 2 

7 2 4- « 2 +2 
cos 0 32 = — 3 25 

2Z 23 s 3 


3 3 = -y3 sec P 
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m i ii mu 


t 

t 


W-, 

y 2 = -it + x 2 tan P 

U 

y3 = " (t + * 3 tan P ) 

32 = y 3 + ^ x 3 " x 2 " y 2 tan P) 


23 


= y| + (x 3 - x 2 + y 3 tan p) 


Z 23 = (x 3 " x 2 )2 + ( y 3 ’ y 2 ^ 1 


cos © 2.4 = cos ® 


41 


14 


z 14 = L 2 + (yj - y 4 )‘ 


2 

2 


cos 0^ = 


L - x 2 
Z 24 


cos ©24 = 


7 2 , 2 2 
^24 + s 2 - t 42 

2Z 2 4 s 2 


t 42 = y| + (L - x 2 - y 2 tan p) 2 


Z 24 = (L - x 2 ) 2 + (y 2 - y 4 )‘ 


W L = Wi + 2L tan p 
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APPENDIX c 


FORTRAN II CODE FOR TUBES AND SLOTS 
by Carl D. Bogart 

The FORTRAN II programs that calculate the vail and exit-plane flux dis- 
tributions and the transmission probabilities follow essentially the same 
format for both the tube and the slot configurations. Thus, vhile both pro- 
grams are included herein, only one set of control vords, one set of problem 
specifications, one flov chart, and one set of sample data are given. 

Control vords: 

IS positive, set first guess equal to YO; zero or negative, read in first 
guess from binary cards 

KO maximum number of iterations on solution 
KPR frequency of intermediate print 

M number of heading cards 

Problem specifications: 


ER 

maximum percentage change for convergence 

H 

reciprocal of mesh size for 

vail 

HH 

reciprocal of mesh size for 

exit plane 

XLOR 

length of tube or slot 


XM 

tan (beta) 


YO 

first guess at solution 
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cc. 


PROGRAM LISTINGS 


C PROGRAM TO CALCULATE SIDEWALL AND EXIT PLANE FLUXFS FOR SLOT 
CHMvnN V,?,w 

DIMFNSI ON Y < 10GO > »W ( 1000 ) *Z ( 1000 > 

1 PF AD INPUT TAPE 7 *101 ,M,KO*KPR , I S 
C M=NUMOFR OF HEADING CARDS 

C KO = MAXIMUM NU^iRFP OF ITERATIONS 

r KPP=FPFOUFNCY OF PRINT 

C IS=0»- CALL FIRST GUESS FRO,' - BINARY CARDS 
C IS = + SFT FIRST Gt)RSS = YO 

DO 2 J = 1*M 

READ INPUT TAPE 7*104 

2 WRITE OUTPUT TAPE 6,104 

READ INPUT TAPE 7 , 1 00 , XM , X LOR , H , ER , YO ,HH 
C X M = T A N ( BFTA } 

C X L 0 P = L / W 
r h=I/M r SH S T 7 p 

r PR = prpcFMTAGF CHANG F FOR C ONVRRGF MCE 
C YQ=F I RST GUESS AT SOLUTION 
C HH = S T EP — S I Z F FOR END 
FL0R=4«*XL0R**2 
F OR = 0 
Knr = n 
HX = 1 ./H 
N=XL0P*H+1 .5 
X s = xm*xm 
C = 1 • + X S 
S C = SO R T F ( C ) 

0SC=1./SC I 

CA=.5*n<;r 1 

TXM=2 • *XM 
FXM=TXM**2 
21 DO 3 J = 1*N 

XJ=J-1 
XJ=XJ/H 

Y( J)=Y0+XJ*.5*( .5-YO) 

K = 1 
x = 0 . 

DO 6 J = 1 » N 
XJ= J-l 
YY=XJ/H 

Z ( J ) = ( 1 • +( X+YY)*TXM+FXM*X*YY )*Y { J ) 

6 Z(J)=Z(J)/(( X-YY ) **2+( 1 .+XM* ( X+YY ) ) **2 ) **1 .5 
5 = 0 * 

ss=o. 

DO 7 J = 2 * M *2 
S=S+Z ( J ) 

7 SS=SS+Z(J+1) 

SOM = ( 7 ( 1 ) -Z ( N >+4.*S + 2.*SS) *. 33353333 /H 
CnN= SC-(XM+X*C)/SQRTF(X*X*C+TXM*X+1« ) 

W(K )=r A*( COM+SOM) 

K = K+1 
X=X+HX 

IF(K-N) 5*5, P 
S CK = 0 • 

DO 10 J= 1 , N 

A= ABSF (W(J)-Y (J) ) / W ( J ) 

IF(CK-A) 9,10,10 

9 CK= A 
K = J 

10 Y ( J ) = W ( J ) 

CK=CK*100. 


r 
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<nc=Knc+i 

KQR = K0R+1 

I F ( KOR-KPR ) 30*31*31 

31 WRITE OUTPUT TAPE 6 , 102 * KOC * K , CK 

WRITE OUTPUT TAPE 6 * 1 0 3 * f Y ( J ) * J= 1 , N ) 

KOR = 0 

30 IF(ER-CK) 18*19,19 

18 I F ( KOC—KO ) 4*19,19 

19 WRITE OUTPUT TAPE 6 , 1 0 2 * KO C > K * CK 
WRITE OUTPUT TAPE 6,105 

WRITE OUTPUT TAPE 6,103, (Y(J) » J = 1 *N ) 

UN=1.+TXM*XL0R 

HN=UN/HH*.5 

NN=HH+1 • 5 

K-l 

X=0. 

39 TX=2.*X 

DO 40 J = 1 * N 
XJ= J-l 
YY=XJ/H 
A= ( XLOR-Y Y ) 

AA = 4 • * A* A 
B= 1 .+TXM*YY 

Z ( J ) =Y ( J) *A* ( ( UN + TX ) / ( AA+ ( B+TX ) **2 ) ** 1 . 5+ ( UN-TX ) / ( AA + (B-TX ) **2 ) 
1 ** 1 • 5 ) 

40 CONTI NUF 
5=0. 

55 = 0 . 

DO 41 J = 2 , N * 2 
S=S+Z ( J ) 

41 5$=$S+Z ( J + l ) 

SOM= ( Z { 1 ) -Z { N ) +4.*S+2.*SS )*.666666666/H 

C 0 N = • 5 * ( ( 1 • -TX ) /SQRTF ( FLOR+ ( 1 • -TX ) **2 ) + ( 1. +TX ) /SQRTE ( FLOR+ ( 1.+ 

1 TX ) **2 ) ) 

W( K) =SOM+cnN 
X= X+HN 
K = K + 1 

IF(K-NN) 39,39,42 

42 WRITE OUTPUT TAPE 6,106 

WRITE OUTPUT TAPE 6 , 1 o 3 , ( W ( J ) , J = 1 , N N ) 

S=0. 

SS= 0 • 

DO 43 J=2*NN*2 
S= 5+W ( J ) 

43 SS=$S+W(J+1) 

SOM= .6666 6666 6* ( W ( 1 ) -W ( NN ) +4. *5+2 . *.SS ) *HN 
WRITE OUTPUT TAPF 6,lO7,S0M 
CALL BCDUMP ( Y ( N ) , Y { 1 ) ) 
on m l 

inn FORMAT (7F1C. 5 ) 

ini FORMAT (1415) 

102 FORMAT ( 2 1HC I TER AT I (IN/ K/ EPS I LON 2 16 ,F8 .3 ) 

103 FORMAT (1H08E15.6 ) 

1D4 FORMAT ( 72FI 

1 ) 

105 FORMAT i 22H0FLUX ON SIDE WALLS ) 

1 06 FORMAT (16H r FLUX OUT F ND ) 

107 FORMAT { 16HOT0TAL OUT FND F10.6) , 

fm n 



C PROGRAM TO CALCULATE SIOFWALL AND FXIT PLANF FLUXES FOR TUBFS 
COMMON Y * Z » W 

D I MENS I ON Y ( ? 0 0 0 ) ,7 ( 2^<'0' ) , W( 2000 ) 

1 RF AD INPUT TAPE 7*101*M»Kn,K n R*I c i 

c M=NUMPFR of HEADING CAPD^ 

C KO=MAXIMUM NUMBER OF ITERATIONS 
C K°R = F REOUFNC Y OF PRINT 

c 1^ = 0, CALL FIRST GUESS F- RUM BINARY CARDS 
C I S = + , S E T FIRST GUFSS = YO 
on 7 J=1»M 

READ INPUT TAPE 7,104 
7 WRITE OUTPUT TAPE 6,104 

PFAD INPUT TAPE 7 , 1 0 0 , XM , X L (1 R , H , FR , YO , HH 
C XM= TAN ( B FT A ) 

C X LOR = L / R 
C H= 1 /M F SH SIZE 

c fr=percfntagf changf for convfrgfnof 

C Y0=FIRST GUESS AT SOLUTION 
C HH = ST C: P — SI Z E FOR END 

KOR = 0 

0=1 .+XM*XM 
SC=SORTF ( C ) 

0SC=1./SC 

<nr = o 

N=XL0R*H+1 . 5 
I F ( IS) 61 ,61 ,6 

61 CALL BCREAD(YfN) , Y { 1 ) > 

GO TO 8 

60 DO 99 J= 1 , N 

99 Y { J ) = Y 0 

P K = 1 

N S = ] 

xx = n. 

9 A=1.+XM*XX 

A A = • 5 / A 
DO 10 J = 1 * N 
X J= J- 1 
X= X J / H 
P=1 .+XM*X 
AL = ABSF ( X-XX ) 

Z ( J ) = AL* { ( AL*AL*C ) +6 .*A*B ) 

I o Z ( J ) = Y ( J)* ( nsc-z ( J ) / ( ( AL*AL*C ) + 4.*A*L3 ) **1 .6 ) 

IF(NS) 11,12,12 

II JJ=XX*H+1.5 

SOM=Z ( 1 ) + • 5*7 {JJ-1 )-Z (JJ)+#5*Z(JJ+1 )-Z ( N ) 

GO TO IS 

1 ? S0M=7 ( 1 ) - 7 ( N ) 

is 5 = ". 

ss=o. 

DO 14 J = 2 , N , 2 
S=S+Z { J ) 

14 SS=SS+7(J+1) 

SDM= ( S0M+4.*S+2.*SS ) / ( S • *H ) 

C 0 N = ( ( C*XX*XX )+3 ,*XM*XX+2 . ) /SORTF ( { C *X X * X X ) +4 . * A ) 

i -( xx+xm*a ) * nsc 

w ( k ) = ( cnw+snM ) *aa 

NS=-NS 

XK = K 

XX= XK / H 

K = K + 1 

IF(K-N) 9,9,15 



15 CK = 0 • 

no i7 J=i »n 

A = A 6 S F ( w ( J ) — y ( J ) )/W{J) 

IF(CK-A) 16*17 >17 

16 CK = A 
K = J 

17 Y(J)=W(J) 

CK = CK*1 00 * 

KOC=KOC+1 

i^op=kor+i 

I F C KOR-KPR ) 30*31 *31 

31 WRITF OUTPUT TAPF 6 * 1 02 * KOC * K * CK 

WR1TF OUTPUT TAPE 6 * 1 03 * ( Y ( J ) * J = 1 ♦ N ) 

K 0 R = 0 

30 IF(ER-CK) 18*19*19 

IP TF(KOC-KO) 8*19*19 

19 WRITE OUTPUT TA D E 6 * 102 * KOC * K * CK 

WRITE OUTPUT TAPE 6*105 
WRITF OUTPUT TAPF 6 * 1 03 * ( Y ( J ) * J = 1 , N ) 

XK = XM*Xt_nR + l . 

NN=HH+1 .5 
HH=XK/HH 
K = 1 
X = 0 • 

11? DO 110 J=1,N 

XJ=J-1 
0=X J/H 
B=1 .+XM*Q 
C = XLOR— 0 

Z ( J ) =Y ( J )*0*C*( (C*C+b*B+X*X ) *XK-2 . *B*X*X ) 

110 Z ( J)=Z ( J)/(X**4+2.*(C*C-B*B)*X*X+(C*C+B*B> **2)**1.5 
SOM = Z C 1 )-Z (N) 

S=0. 

?s=0 # 

no 111 J = 2 * N * 2 
5= S + Z ( J) 

111 S S = SS + Z ( J+ 1 ) 

SOM = ( S0M+4.*S+2.*SS)*. 666666666 /H 

C0N = • 5* ( 1 • - ( ( XLOR **2—1 . ) + X*X ) /SQRTF { X**4+2 .* ( XL0P**2-3 • ) *X*X + 
1 (XL0R**2+1.)**2) ) 

W(K)=SQM+cOM 
X= X + HH 
K=^ + l 

IF (K-NN) 112*112*113 

113 W p I T F OUTPUT TAPF 6*106 

WRITF OUTPUT TAPE 6 * 1 0 3 * ( W ( J ) * J = 1 * NN ) 

X = 0 # 

DO 1 1 A J = 1 * N N 
Z( J)=W( J)*X 

114 X= X+HH 
o = o # 

ss= o • 

DO 11S J= ? * NN * 2 
S= S+Z ( J ) 

115 SS=SS+Z(J+1) 

S 0 M = H H * ( Z ( 1 ) -Z ( NN ) +4, *S + 2 .*SS ) *.6666 666 66 
WRITE OUTPUT TAPF 6 * 107 * r > 0 M 
CALI BC DUMP ( Y t N ) * Y f 1 ) ) 
an to i 

10" FORMAT (7F10.5 ) 

1"! FORMAT ( ]4I5 ) 
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1"2 

103 

1^4 

1 

106 

107 


*DATA 

1 

XM = 
• IB 


F flRM AT ( 2 1H0 I TERAT I ON/K /FPS I LON 2 I 6 *K8*3 ) 
FORMAT ( 1H08F15.6 ) 

FflPMAT ( ^ 2 H 

) 

FORMAT (22HOFLUX ON M DF WALLS ) 

FORMAT ( 16 H 0 FLUX OlJT F NO ) 

FORMAT ( 16H0T0TAL OUT END F10.6) 

FNO 


40 41 1 

18*XLQR=2.*H=20.fER=.01*YO=0.»HH=?00. 

2.p 20. .01 0. 200. 
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TABLE I. - TOTAL TRANSMISSION PROBABILITI 


(a) Tubes 


Length to 
inlet -radius 
ratio, 
L/*l 

Wall half -angle, (3, deg 

Diverging walls 

Parallel 

walls 

Converging walls 

75 

60 

45 

30 

20 

10 

5 

0 

-5 

-10 

-20 

-30 

-45 

-60 

Total transmission probability, P^. 

0.5 

1 

2 

4 

8 

16 

0. 999 
.999 
.999 
.999 

0.996 

.996 

.996 

.995 

0.984 0.953 0.917 

. Q7Q . . R7 A 

0. 869 
.792 
.708 
.638 
.594 
.573 

0. 834 
.730 
.608 
.495 
.410 
.359 

0. 801 
.671 
.513 
.355 
.223 
.130 

1 

0. 763 
.606 
.408 
.204 
.044 

0.712 0.609 
.516 .347 

<?7n nco 

0. 475 
.161 

0.245 

0. 022 

.. Q7£ . . Q1 ft 

.834 

.807 

.795 

.791 



.975 

.975 

.910 

.908 

.907 

.045 


— 





.... 










(b) Slots 


Length to 
inlet -width 
ratio, 
L/Wi 


Wall half -angle, (3, deg 


Diverging walls 

Parallel 

walls 

Converging walls 

75 

60 

45 

30 

20 

10 

5 

0 

-5 

-10 

1 

to 

o 

-30 

-45 

-60 


Total transmission probability, Bfc 


0. 25 

0.999 

0. 9975 

0.990 0.973 

0.953 

0. 927 

0. 908 

0.889 

0. 868 

0.838 

0. 775 

0. 686 

0. 493 

0.149 

.5 

.999 

.9970 

.986 

.958 

.924 

.876 

.840 

.804 

.763 

.704 

.578 

.394 



1 

.999 

» 9966 

.983 

.944 

.891 

.812 

.748 

.684 

.607 

.493 

.237 





2 

.999 

.9964 

.981 

.935 

.865 

.750 

.648 

.541 

.405 

.188 





4 

.999 

.996 

.981 

.930 

.848 

.701 

,560 

.398 

.169 







8 



.981 

.928 

.839 

.671 

.495 

.271 

t 




i 




CO 



(A 

O 


TABLE II. - DIRECT TRANSMISSION PROBABILITY^ 


(a) Tubes 


Length to 
inlet-radius 
ratio, 
L/Ri 






Wall half- angle, f 

S, deg 






Diverging walls 

Parallel 

walls 

Converging walls 




75 

60 

45 

30 

20 

10 

5 

0 

-5 

-10 

-20 

-30 

-45 

-60 






Direct transmission probability, ] 

?d 





0.5 

1 

2 

A 

0.966 
.955 
. 946 
. 940 

0.910 
.865 
.823 
. 790 

0.849 
. 764 
.675 
.604 
.557 

0. 780 
.649 
.510 
.399 
.330 
.292 

0. 727 
.563 
.391 
.262 
.188 
.151 

0.673 

.478 

.282 

.149 

.0840 

.0553 

0.639 

.426 

.220 

.0934 

.0398 

.0198 

0.610 
.382 
.172 
.0557 
.0152 
. 00388 

0.578 

.337 

.126 

.0265 

.0020 

0.538 

.281 

.0768 

.0046 

1 

0.460 

.184 

.0155 

1 

0. 362 
.0843 

0. 191 

0. 0179 

8 

16 













(b) Slots 


Length to 
inlet-width 
ratio, 

l/w-l 

Wall half- angle, p, deg 

Diverging walls 

Parallel 

walls 

Converging walls 


75 

60 

45 

30 

^0 

10 

5 

0 

-5 

-10 

-20 

-30 

-45 

-60 

Direct transmission probability, 

0. 25 

.5 

1 

2 

4 

8 

0. 983 
.977 
. 972 
.969 

0.954 

.930 

.907 

.889 

0.921 
.874 
.822 
. 777 
.746 

0.883 

.805 

.714 

.632 

.574 

.540 

0. 853 
.750 
.625 
.512 
.434 
.389 

0. 820 
.691 
.531 
. 386 
.290 
.235 

0. 799 
.652 
.469 
.306 
.199 
.141 

0.781 

.618 

.414 

.236 

.123 

.0623 

0. 761 
.581 
.356 
.163 
.0446 

0. 733 
.530 
.277 
.0678 

0.678 

.429 

.124 

0.602 

.290 

0.437 

0.134 


















TABLE III 


WALL FLUX DISTRIBUTIONS 


Length to Axial 
inlet-radius dis- 
ratio, tance, 

L/R-l x/L 


(a) Tubes 

Wall half-angle, £, deg 

Parallel 
walls 

20 | 10 | 5 0 -5 

Flux ratio, n 2 /n]_ 

0.0176 0.0727 0.161 0.286 0.388 0.491 0.553 

.0124 .0614 .145 .265 .367 .470 .532 

.00911 .0523 .130 .247 .346 .449 .511 

.00690 .0449 .118 .230 .327 .428 .490 

.00535 .0388 .106 .214 .308 .408 .470 

.00424 .0338 .0964 .197 .290 .388 .449 

.00341 .0296 .0876 .185 .273 .369 .429 

.00279 .0260 .0798 .172 .257 .350 .409 

.00231 .0230 .0727 .160 .242 .332 .389 

.00194 .0204 .0664 .149 .227 .314 .370 

.00164 .0182 .0608 .139 .213 .297 .350 

0.0176 0.0734 0.166 0.302 0.420 0.541 0.614 0.674 0.732 

.00913 .0530 .135 .264 .380 .502 .578 .640 .700 

.00537 .0396 .111 .232 .344 .465 .542 .605 .668 

.00343 .0303 .0927 .204 .311 .430 .506 .570 .635 

.00233 .0237 .0779 .179 .280 .396 .471 .535 .601 

.00166 .0189 .0659 .158 .253 .363 .437 .500 .566 

.00122 .0153 .0562 .140 .228 .332 .403 .465 .529 

.000924 .0126 .0482 .124 .205 .303 .371 .430 .492 

.000717 .0104 .0415 .109 .184 .275 .339 .395 .454 

.000568 .0088 .0360 .0969 .165 .249 .308 .360 .416 

.000457 .00741 .0313 .0859 .148 .225 .278 .326 .377 

0.0176 0.0737 0.168 0.315 0.450 0.595 0.685 0.758 0.827 

.00537 .0399 .114 .247 .378 .529 .626 .707 .786 

.00233 .0241 .0811 .196 .318 .468 .569 .656 .744 

.00122 .0157 .0595 .157 .269 .412 .514 .604 .698 

.000722 .0108 .0449 .127 .227 .362 .461 .552 .650 

.000462 .00775 .0347 .104 .192 .316 .411 .500 .598 

.000314 .00576 .0273 .0854 .163 .275 .363 .448 .544 

.000223 .00440 .0218 .0708 .138 .237 .317 .396 .487 

.000164 .00343 .0177 .0590 .116 .203 .274 .344 .427 

.000124 .00273 .0145 .0492 .0982 .172 .233 .293 .364 

. 0000965 . 00220 . 012 . 0412 . 0824 . 144 . 194 . 242 . 298 

0.0176 0.0737 0.169 0.322 0.469 0.638 0.748 0.836 0.917 

.00234 .0242 .0824 .204 .344 .526 .656 .769 .879 

.000723 .0109 .0463 .137 .257 .433 .572 .701 .834 

.000315 .00587 .0288 .0959 .195 .357 .495 . 633 .786 

.000165 .00354 .0192 .0698 .151 .295 .426 .566 .733 

.0000976 .00231 .0134 .0522 -118 .242 .363 .500 .674 

.0000625 .00160 .00980 .0399 .0935 .198 .306 .433 .607 

.0000425 .00115 .00734 .0309 .0740 .161 .253 .366 .531 

.0000302 .000856 .00564 .0242 .0587 .129 .204 .298 .443 

.0000222 .000654 .00441 .0192 .0464 .101 .159 .230 .341 

.0000168 .000511 .00349 .0152 .0364 .0770 .117 .163 .224 

0.170 0.325 0.478 0.665 0.793 0.899 0.983 

.047 .140 .270 .476 .650 .816 .966 

.020 .074 .168 .349 .533 .734 .945 

.010 .044 .112 .262 .437 .655 .920 

.0061 .029 .078 .200 .357 .577 .889 

.004 .020 .056 .153 .290 .500 .849 

.0027 .014 .042 .118 .231 .422 .795 

.0019 .010 .031 .090 .181 .344 .721 

.0014 .0079 .023 .067 .136 .265 .615 

.0011 .0060 .017 .049 .0959 .184 .452 

.00085 .0045 .0129 .034 .0608 .101 .173 

0.326 0.482 0.678 0.820 0.942 

.0750 .173 .375 .594 .847 

.0302 .0849 .232 .442 .757 

.0157 .0491 .155 .335 .670 

.00936 .0312 .108 .256 .584 

.00606 .0210 .0769 .195 .499 

.00415 .0147 .0554 .147 .413 

.00294 .0105 .0399 .108 .327 

.00213 .00753 .0282 .0764 .241 

.00157 .00540 .0192 .0494 .152 

.00117 .00384 .0122 .0267 .0579 


0.604 0.653 

.583 .633 

.562 .613 

.542 .593 

.521 .572 

.500 .551 

.479 .530 

.458 .509 

.438 .487 

.417 .465 

.396 .443 


Diverging walls 
75 I 60 I 45 I 30 | 
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TABLE III 


Concluded. WALL FLUX DISTRIBUTIONS 


(b) Slots 


Length to 

Axial 





Wall half-angle. 

P, deg 

Inlet -width 

dis- 









ratio. 

tance. 



Diverging walls 



Parallel 

L/Wi 

x/L 








walls 



75 

1 

60 

| 45 

| 30 

1 20 

I 10 

1 5 

0 








Flux r 

atio, n 

2/ n l 

0.25 

0 

0.0175 

0.0716 

0.156 

0.271 

0.365 

0.458 

0.515 

0.561 


.1 

.0147 

.0656 

.147 

.261 

.353 

.446 

.503 

.549 


.2 

.0124 

.0602 

,139 

.250 

.342 

.434 

.491 

.537 


.3 

.0107 

.0555 

.132 

.240 

.331 

.422 

.478 

.524 


.4 

.00927 

.0513 

.125 

.231 

.320 

.411 

.466 

.512 


.5 

.00812 

.0475 

.118 

.222 

.309 

.399 

.454 

.500 


.6 

.00717 

.0441 

.112 

.213 

.299 

.388 

.442 

.488 


.7 

.00638 

.0411 

.106 

.205 

.289 

.377 

.431 

.476 


.8 

.00571 

.0383 

.101 

.197 

.279 

.366 

.419 

.463 


.9 

.00514 

.0358 

.0961 

.189 

.270 

.355 

.407 

.451 


1.0 

.00465 

.0336 

.0915 

.182 

.261 

.344 

.396 

.439 

0.5 

0 

0.0176 

0.0724 

0,160 

0.285 

0.389 

0.495 

0.560 

0.613 


.1 

.0125 

.0611 

.144 

.264 

.367 

.473 

.538 

.591 


.2 

.00931 

.0522 

.129 

.245 

.345 

.450 

.515 

.569 


.3 

.00721 

.0450 

.117 

.227 

.325 

.428 

.492 

.546 


.4 

.00574 

.0392 

.106 

.211 

.305 

.406 

.470 

.523 


.5 

.00469 

.0344 

.0960 

.196 

.287 

.385 

.447 

.500 


.6 

.00390 

.0304 

.0875 

.182 

.269 

.364 

.425 

.477 


.7 

.00329 

.0271 

.0800 

.170 

.253 

.345 

.404 

.454 


.8 

.00281 

.0243 

.0734 

.158 

.238 

.326 

.382 

.431 


.9 

.00246 • 

.0218 

.0675 

.147 

.223 

.307 

.362 

.409 


1.0 

.00213 

.0198 

.0622 

.138 

.210 

.290 

.342 

.387 

1 

0 

0.0176 

0.0730 

0.164 

0.300 

0.418 

0.544 

0.623 

0.688 


.1 

.00932 

.0528 

.133 

.261 

.377 

.503 

.584 

.653 


.2 

.00576 

.0398 

.110 

.227 

.338 

.464 

.546 

.616 


.3 

.00391 

.0310 

.0917 

.199 

.304 

.425 

.507 

.578 


.4 

.00283 

.0249 

.0776 

.175 

.273 

.389 

.468 

.539 


.5 

.00214 

.0204 

.0664 

.154 

.245 

.355 

.431 

.500 


.6 

.00168 

.0170 

.0573 

.137 

.220 

.323 

.395 

.461 


.7 

.00135 

.0144 

.0500 

.122 

.198 

.293 

.360 

.422 


.8 

.00111 

.0123 

.0440 

.109 

.178 

.265 

.327 

.384 


.9 

.000934 

.0107 

.0389 

.0975 

.161 

.240 

.296 

.347 


1.0 

.000794 

.00937 

.0347 

.0878 

.146 

.217 

.267 

.312 

2 

0 

1 0.0176 

0.0732 

0.166 

0.310 

0.442 

0.590 

0.688 

0.770 


.1 

.00577 

.0401 

.112 

.239 

.367 

.520 

.625 

.719 


.2 

.00284 

.0252 

.0803 

.188 

.305 

.455 

.564 

.665 


.3 

.00169 

.0173 

.0602 

.151 

.255 

.397 

.505 

.610 


.4 

.0011? 

.0127 

.0468 

.123 

.215 

.345 

.450 

.555 


.5 

.000800 

.00968 

.0374 

.102 

.183 

.300 

.398 

.500 


.6 

.000600 

.00765 

.0306 

.0858 

.156 

.261 

.349 

.445 


.7 

.000467 

.00621 

.0255 

.0729 

.134 

.225 

.304 

.390 


.8 

.000374 

.00514 

.0216 

.0625 

.115 

.194 

.261 

.335 


1 .9 

.000306 

.00433 

.0185 

.0541 

.0996 

.166 

.222 

.281 


1.0 

.000256 

.00370 

.0161 

.0471 

.0865 

.142 

.186 

.230 

4 

0 

0.0176 

0.0733 

0.167 

0.315 

0.457 

0.624 

0.740 

0.843 


.1 

.00284 

.0253 

.0815 

.195 

.325 

.502 

.640 

.773 


.2 

.00112 

.0128 

.0482 

.132 

.240 

.405 

.549 

.701 


.3 

.000602 

.00779 

.0321 

.0950 

.183 

.331 

.471 

.632 


.4 

.000376 

.00528 

.0231 

.0720 

.144 

.273 

.403 

.565 


.5 

.000258 

.00383 

.0175 

.0565 

.116 

.226 

.343 

.500 


.6 

.000188 

.00292 

.0137 

.0454 

.0946 

.187 

.289 

.433 


.7 

.000143 

.00230 

.0111 

.0372 

.0778 

.154 

.240 

.366 


.8 

.000113 

.00187 

.00914 

.0309 

.0645 

.127 

.195 

.298 


.9 

.0000914 

.00154 

.00767 

.0260 

.0538 

.103 

.154 

.226 


1.0 

. 0000755 

.00130 

.00653 

.0221 

.0451 

.0829 

.117 

.156 

8 

0 



0.167 

0.317 

0.464 

0.644 

0.776 

0.899 


.1 



.0486 

.135 

.251 

.439 

.613 

.810 


.2 



.0236 

.0762 

.159 

.317 

.492 

.725 


.3 



.0143 

.0501 

.111 

.240 

.400 

.647 


.4 



.00976 

.0358 

.0824 

.187 

.327 

.571 


.5 



.00714 

.0270 

.0634 

.148 

.267 

.496 


.6 



.00547 

.0211 

.0498 

.117 

.216 

.422 


.7 



.00433 

.0168 

.0397 

.0932 

.172 

.347 


.8 



.00352 

.0137 

.0320 

.0733 

.133 

.269 


.9 



.00292 

.0113 

.0261 

.0569 

.0977 

.187 


1.0 



.00246 

.00949 

.0213 

.0436 

.0673 

.0990 


Converging walls 


-5 -10 -20 -30 -45 -60 


0.607 

.595 

.583 

.571 

.558 

.546 

.533 

.521 

.508 

.496 

.483 


0.662 

.650 

.639 

.627 

.615 

.602 

.590 

.577 

.564 

.551 

.538 


0.752 

.741 

.731 

.719 

.708 

.696 

.683 

.670 

.657 

.643 

.629 


0.838 

.830 

.820 

.811 

.800 

.789 

.777 

.765 

.751 

.737 

.722 


0.938 

.933 

.927 

.920 

.913 

.904 

.894 

.883 

.871 

.856 

.840 


0.995 

.994 

.993 

.992 

.990 

.988 

.985 

.979 

.971 

.956 

.927 


0.666 

.645 

.623 

.601 

.578 

.554 

.530 

.506 

.481 

.457 

.433 


0.730 

.710 

.690 

.668 

.645 

.622 

.597 

.571 

.545 

.518 

.491 


0.831 

.816 

.799 

,780 

.760 

.737 

.713 

.686 

.657 

.625 

.592 


0.924 

.915 

.904 

.891 

.876 

.858 

.837 

.811 

.780 

.743 

.697 


0.754 

.723 

.689 

.653 

.615 

.575 

.534 

.491 

.448 

.404 

.360 


0.832 

.808 

.780 

.750 

.715 

.677 

.635 

.589 

.538 

.483 

.425 


0.949 

.940 

.928 

.914 

.896 

.874 

.844 

.804 

.748 

.665 

.545 


0.854 

.817 

.775 

.730 

.682 

.629 

.572 

.509 

.439 

.362 

.281 


0.951 

.937 

.919 

.897 

.871 

.839 — 

.798 

.743 

.666 

.547 

.358 


0.949 

.923 

.892 

.'857 

.817 

.770 

.713 

.641 

.547 

.413 

.212 
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Length to 
[inlet-radius 
ratio, 

L/R;l 


Radial 

dis- 

tance. 


0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 


.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 

0 

.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.95 

1.0 


TABLE IV. - EXIT- PLANE FLUX DISTRIBUTIONS 
(a) Tubes 

Wall half-angle, (3, deg 


0.802 

.780 

.699 

.515 

.267 

.113 

.0513 

.0265 

.0152 

.00946 

.00763 

.00461 

0.504 

.450 

.303 

.153 

.0691 

.0329 

.0171 

.00971 

.00530 

.00378 

.00308 

.00209 

0.204 

.164 

.0912 

.0426 

.0200 

.0100 

.00546 

.00318 

.00196 

.00127 

.00104 

.000762 

0.0609 
.0471 
.0251 
.0118 
.00564 
.00269 
.00159 
.000333 
.000579 
.000377 
. 000308 
. 00023Q 


0.810 

.802 

.773 

.717 

.619 

.474 

.315 

.193 

.119 

.0758 

.0615 

.0331 

0.520 

.502 

.450 

.368 

.274 

.188 

.125 

.0823 

.0552 

.0377 

.0313 

.0191 

0.220 

.208 

.176 

.135 

.0973 

.0673 

.0459 

.0314 

.0217 

.0152 

.0128 

.00867 

0.0694 

.0653 

.0548 

.0420 

.0304 

.0213 

.0148 

.0103 

.00719 

.00510 

.00431 

.00318 


Diverging walls 
45 1 30 I 20 1 

i° [ 

5 

Parallel 

walls 

0 

.824 

0.846 

0.865 

Flux rat 
0.884 

io, n4/ 
0.896 

n l 

0.906 

.819 

.843 

.862 

.882 

.894 

.904 

.803 

.832 

.854 

.876 

.890 

.901 

.772 

.813 

.840 

.866 

.881 

.894 

.720 

.782 

.819 

.851 

.869 

.883 

.641 

.735 

.786 

.828 

.850 

.868 

.531 

.669 

.740 

.795 

.824 

.846 

.406 

.582 

.676 

.750 

.788 

.816 

.293 

.481 

.597 

.691 

.740 

.776 

.206 

.380 

.506 

.617 

.676 

.721 

.172 

.333 

.458 

.473 

.637 

.686 

.0903 

.180 

.259 

.343 

.395 

.439 

.550 

0.600 

0.648 

0.700 

0.734 

0.762 

.541 

.595 

.644 

.698 

.732 

.761 

.513 

.579 

.634 

.691 

.726 

.756 

.468 

.553 

.616 

.678 

.717 

.749 

.408 

.517 

.590 

.661 

.703 

.737 

.340 

.471 

.558 

.638 

.685 

.723 

.271 

.419 

.519 

.610 

.662 

.704 

.210 

.363 

.473 

.575 

.634 

.680 

.159 

.307 

.423 

.534 

.599 

.650 

.119 

.253 

.369 

.486 

.555 

.611 

.102 

.226 

.339 

.456 

.528 

,586 

.0599 

.133 

.205 

.286 

.340 

.386 

.252 

0.313 

0.380 

0.463 

0.522 

0.576 

.246 

.310 

.378 

.462 

.521 

.575 

.230 

.301 

.371 

.457 

.518 

.572 

.205 

.286 i 

.360 

.449 

.511 

.566 

.176 

.266 

.346 

.439 

.502 

.559 

.146 

.243 

.328 

.425 

.491 

.549 

.118 

.217 

.306 

.408 

.476 

.536 

.0930 

.191 

.282 

.387 

.458 

.520 

.0724 

.164 

.255 

.362 

.435 

.499 

.0556 

.138 

.226 

.333 

.407 

.472 

.0484 

.125 

.209 

.315 

.389 

.455 

.0314 

.0797 

.134 

.206 

.258 

.307 

).0879 

0.128 

0.179 

0.257 

0.325 

0.396 

.0859 : 

.126 

.178 

.256 

.325 

.396 

.0802 

.123 

.175 1 

.254 

.322 

.394 

.0717 

.117 

.170 

.250 

.319 

.390 

.0619 

.109 

.164 

.245 

.314 

.385 

.0518 

.101 

.156 

.238 : 

.307 

.379 

.0423 

.0308 

.147 

.229 

.298 

.371 

.0339 

.0805 

.136 

.218 

.288 

.360 

.0268 

.0700 

.124 

.205 

.275 

.347 

.0203 

.0596 

.111 

.190 

.258 

.329 

.0183 

.0544 

.103 

.180 

.247 

.318 

.0130 

.0373 

.0701 

.122 

.169 

.219 

). 0264 

0.0431 

0. 0684 

0.118 

0.174 

0.249 

.0258 

.0427 

.0681 

.118 

.173 

.249 

.0241 

.0416 

.0671 

.117 

.172 

.248 

.0217 

.0397 

.0654 

.115 

.170 

.246 

.0188 

.0373 

.0631 

.112 

.168 

.243 

.0159 

.0344 

.0602 

.109 

.164 

.239 

.0131 

.0313 

.0568 

.106 

.160 

.234 

.0106 

.0279 

.0529 

.101 

.154 

.227 

.00845 

.0245 

.0486 

.0953 

.148 

.219 

.00667 

.0211 

.0437 

.0885 

.139 

.209 

.00590 

.0194 

.0411 

.0845 

.134 

.200 

.00446 

.0141 

.0292 

.0590 

.0931 

.140 


0.0127 

0.0221 

0.0446 

0.0782 

0.145 


.0126 

.0220 

.0445 

.0781 

.145 


.0123 

.0217 

.0441 

.0776 

.144 


.0118 

.0212 

.0435 

.0768 

.143 


.0111 

.0205 

.0426 

.0757 

.141 


.0103 

.0196 

.0415 

.0742 

.139 


.00936 

.0185 

.0401 

.0723 

.136 


.00840 

.0173 

.0384 

.0699 

.132 


.00743 

.0160 

.0364 

.0670 

.128 


.00645 

.0145 

.0339 

.0632 

.122 


.00596 

.0137 

.0325 

.0609 

.117 


.00457 

.0101 

.0233 

.0430 

.0824 


Cunv^rging walls 


-10 -20 -30 -45 -60 


0.916 

.914 

.911 

.906 

.897 

.ee4 

.867 

.842 

.809 

.763 

.732 

.483 

0.791 

.790 

.786 

.780 

.771 

.759 

.744 

.724 

.699 

.665 

.642 

.434 

0.633 

.632 

.629 

.625 

.619 

.610 

.599 

.585 

.567 

.542 

.525 

.361 

0.484 

.483 

.481 

.478 

.474 

.467 

.459 

.449 

.436 

.418 

.405 

.282 

0.375 

.375 

.373 

.371 

.367 

.363 

.357 

.349 

.339 

.316 

.283 

.219 


0.927 

.927 

.924 

.920 

.913 

.903 

.890 

.871 

.846 

.809 

.784 

.539 

0.826 

.825 

.823 

.818 

.812 

.803 

.791 

.776 

.756 

.728 

.709 

.495 

0.706 

.705 

.703 

.700 

.695 

.688 

.679 

.667 

.652 

.630 

.614 

.434 

0.613 

.613 

.611 

.608 

.604 

.598 

.591 

.581 

.568 

.546 

.509 

.379 


0.947 
.946! 
.945 i 
.942! 
.938 
.932 
.924 
.913 
.898 
.875 
.859 
.631 

0.883 

.883 

.881 

.878 

.875 

.870 

.863 

.854 

.841 

.823 

.809 

.600 

0.827 

.827 

.826 

.823 

.820 

.816 

.810 

.802 

.791 

.774 

.756 

.566 


.965 

.965 

.964 

.963 

.961 

.958 

.954 

.948 

.940 

.928 

.918 

.725 

0.934 

.934 

.933 

.932 

.930 

.928 

.924 

.919 

.912 

.902 

.894 

.708 


.986 

.965 

.965 

.985 

.984 

.983 

.982 

.980 

.978 

.974 

.971 

.842 


|0.997 
.397 
.997 
.936 
.996 
.996 
.996 : 
,996 
.995 
.994 | 
.991 
.928 
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TABLE IV 
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EXIT- PLANE FLUX DISTRIBUTIONS 
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I 
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1 
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[ 
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5 

0 

-5 
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Flux 
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0.25 

0 

0.896 

0.901 

0.909 

0.920 

0.929 

0.939 

0.944 

0.949 

0.954 

0.960 

0.969 

0.978 

0.990 

0.997 


.1 

.877 

.893 

.904 

.917 

.927 

.937 

.943 

.948 

.953 

.959 

.969 

.978 

.990 

.997 


.2 

.805 

.870 

.890 

.908 

.921 

.932 

.939 

.945 

.950 

.957 

.967 

.978 

.989 

.997 


.3 

.627 

.820 

.864 

.892 

.909 

.924 

.932 

.939 

.945 

.953 

.965 

.976 

.989 

.997 


.4 

.361 

.728 

.818 

.865 

.890 

.910 

.921 

.930 

.938 

.947 

.961 

.974 

.988 

.997 


.5 

.175 

.584 

.744 

.823 

.861 

.890 

.905 

.916 

.927 

.938 

.956 

.972 

.988 

.997 


.6 

.0918 

.414 

.636 

.761 

.818 

.860 

.881 

.897 

.911 

.926 

.949 

.968 

.987 

.997 


.7 

.0541 

.274 

.507 

.676 

.759 

.819 

.848 

.870 

.889 

.910 

.939 

.963 

.985 

.996 


.8 

.0351 

.183 

.384 

.575 

.682 

.763 

.803 

.832 

.858 

.886 

.925 

.956 

.983 

.996 


.9 
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.128 

.288 

.473 

.594 

.693 

.745 

.783 

.817 

.854 
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.946 

.980 

.993 
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.109 

.250 

.425 

.547 
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.710 

.752 

.791 

.833 
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.939 

.978 
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.453 
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.582 
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.747 

.851 

.930 
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0 

0.710 

0.722 

0.743 

0.774 
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0.828 

0.845 
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0.892 
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.1 

.652 

.704 

.733 

.768 

.797 

.826 

.844 

.858 

.873 

.892 

.922 

.951 
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.480 

.647 

.705 

.752 

.786 

.819 

.838 

.854 

.870 

.889 

.920 

.950 

— 
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.279 

.555 

.656 

.725 

.768 

.807 

.829 

.847 

.864 

.885 

.918 

.949 




.4 

.149 

.439 

.589 

.687 

.743 

.790 

.816 

.836 

.856 

.879 

.915 

.948 




.5 

.0835 

.326 

.509 

.638 

.709 

.767 

.798 

.822 

.845 

.871 

.910 

.945 




.6 

.0505 

.235 

.425 

.581 

.668 

.739 

.776 

.804 

.830 

.860 

.904 

.943 




.7 

.0327 

.170 

.347 

.518 

.621 

.705 

.749 

.782 

.812 

.847 

.897 

.939 




.8 

.0225 

.125 

.279 

.453 

.568 

.664 

.716 

.754 

.790 

.829 

.887 

.934 




.9 

.0161 

.9833 

.223 

.390 

.510 

.617 

.675 

.719 

.761 

.807 

.873 

.926 

— 
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.0813 
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.479 
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.651 

.698 
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1.0 
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.0526 

.124 

.223 

.304 

.385 

.435 

.475 

.516 

.567 

.651 

.737 



1 

0 

0.452 

0.470 

0.502 

0.554 

0.603 

0.657 

0.692 

0.723 

0.755 

0.795 

0.864 





.1 

.384 

.450 

.493 

.550 

.600 

.655 

.691 

.722 

.754 

.794 

.864 


— 



.2 

.245 

.397 

.467 

.536 

.592 

.650 

.687 

.719 

.751 

.793 

.863 





.3 

.137 

.325 

.428 

.515 

.578 

.641 

.680 

.714 

.747 

.790 

.861 





.4 

.0770 

.252 

.380 

.487 

.559 

.629 

.671 

.706 

.742 

.785 

.859 





.5 

.0457 

.190 

.329 

.453 

.536 

.613 

.659 

.696 

.734 

.780 

.856 





.6 

.0289 

.142 

.279 

.416 

.509 

.594 

.644 

.684 

.724 

.772 

.851 





.7 

.0192 

.106 

.233 

.377 

.478 

.571 

.625 

.669 

.711 

.763 

.845 





.8 

.0134 
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.193 

.336 

.443 

.544 

.603 
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.695 

.750 

.837 





.9 

.00969 
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.158 
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.512 
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.625 
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.733 
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. 00834 
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.384 
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.188 

.263 

.341 

.391 

.432 

.475 

.529 

.624 




2 

0 

0.247 

0.264 

0.295 

0.349 

0.405 

0.473 

0.522 

0.568 

0.620 

0.694 






.1 

.204 

.252 

.290 

.347 

.403 

.472 

.521 

.567 

.619 

.693 






.2 

.127 

.220 

.275 

.339 

.398 

.469 

.519 

.565 
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.692 



— 



.3 

.0715 

.180 

.252 

.327 

.390 

.463 

.514 

.562 

.615 

.690 






.4 

.0411 
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.226 

.311 

.379 

.456 

.509 

.557 

.611 

.686 






.5 

.0248 

.108 

.197 
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.365 

.446 

.501 

.550 
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.682 






.6 

.0159 
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.169 

.270 

.349 

.434 

.491 

.542 

.598 

.676 






.7 
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.0623 

.143 

.247 

.330 

.420 

.479 

.532 

.589 

.669 






.8 

.00746 

.0478 

.120 

.223 

.309 

.403 

.464 

.519 

.578 

.659 






.9 

.00541 

.0370 

.100 

.199 

.286 

.382 

.446 

.502 

.563 

.646 






.95 

.00466 

.0327 

.0913 
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.273 

.370 
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0.161 

0.201 
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.373 
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(a) Tube. 



Figure 5. - Comparison of transmission probabilities through cylindrical tubes and 
parallel-walled slots. 


59 


I 




Wall half-angle, 

JB. 

deg 



_Wall half-angle, J 

A 

deg 



-Wall half-angle, 


Axial distance, x/L 

(a) Length to inlet-radius ratio, 0.5. (b) Length to inlet-radius ratio, 2. (c) Length to inlet-radius ratio, 16. 

Figure 6. - Variation of flux along wall of convergent and divergent tubes. (Diverging walls; positive converging walls, negative / 3 .) 






(a) Length to inlet-radius ratio, 0.5. 

Figure 9. - Variation of flux across exit 
negative /3.) 


Radial distance, r/R|_ 

(b) Length to inlet-radius ratio, 2. (c) Length to inlet-radius ratio, 16. 

e of convergent and divergent tubes. (Diverging walls, positive /3; converging walls, 
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